Running faster together: huge speed up of thermal ratchets due to hydrodynamic 

coupling 
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We present simulations that reveal a surprisingly large effect of hydrodynamic coupling on the 
speed of thermal ratchet motors. The model that we use considers particles performing thermal 
ratchet motion in a hydrodynamic solvent. Using particle-based, mesoscopic simulations that main- 
tain local momentum conservation, we analyze quantitatively how the coupling to the surrounding 
fluid affects ratchet motion. We find that coupling can increase the mean velocity of the moving 
particles by almost two orders of magnitude, precisely because ratchet motion has both a diffusive 
and a deterministic component. The resulting coupling also leads to the formation of aggregates 
at longer times. The correlated motion that we describe increases the efficiency of motor-delivered 
cargo transport and we speculate that the mechanism that we have uncovered may play a key role 
in speeding up molecular motor-driven intracellular transport. 

PACS numbers: 67.40.Hf, 87.16.Nn, 05.40.-a, 87.16.Wd 



The motion of particles in a fluid is affected by the 
hydrodynamic interaction mediated by the embedding 
medium [![. The dynamic coupling of small suspended 
particles affects, for instance, the average velocity of par- 
ticles driven by a constant force [l[ and the correlation 
spectrum of pairs of freely diffusing particles 0, Q • All 
existing studies show that hydrodynamic coupling can 
increase the speed at which particles move, but the re- 
sulting speed-ups are typically quite modest [3, Q . 

In this Letter we show that hydrodynamic interactions 
can cause a very large speedup of particles that move 
asynchronously by a thermal ratchet mechanism: hydro- 
dynamic coupling increases the speed of such motors by 
up to two orders of magnitude compared to the velocity 
of isolated particles. Physical realizations of such ratchet 
motors, to which we will refer generically as steppers, can 
be created in colloidal systems and may be found in 
molecular motors Q that move along polar biofilaments, 
such as microtubules or actin. Hence, the effect of hy- 
drodynamic coupling on stepping particles is likely to 
be relevant for the understanding of the physical mech- 
anisms underlying intracellular transport processes such 
as cytoplasmic streaming axonal transport [siGlJ and 
membrane-embedded cargo pulling 
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In order to study the behavior of many steppers mov- 
ing along the same filament, we employ a simple model 
that accounts both for the essential features of steppers 
and for the time-dependent hydrodynamics of the embed- 
ding fluid. The moving particles are described using the 
two-state ratchet model Q, a standard, simplified model 
that accounts for the mechanochcmical coupling under- 
lying molecular motor mechanics. In this ratchet model, 
the stepper can be in two different internal states: in state 
1 particles displace under the action of a potential, V{x), 
of period I, which depends on the position, x , along the 
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FIG. 1: Typical trajectory generated by a two-state ratchet 
model for a single stepper. Inset: free-energy landscape. 
In the 7 region of state 1 (power stroke phase) the particle 
experiences the potential Vi. In region e the particle can jump 
to the state 2 with a rate pi2- In state 2 (diffusive phase) the 
particle diffuses and can revert to state 1 with a rate P2i- 
Parameter values used in the simulations: 7 = 0.75, e = 0.01, 
S = 0.2, AVi = 200 (maximum energy difference); lengths in 
units of the ratchet period I and energy in units of ksT. 



filament axis and that we consider piecewise-linear for 
simplicity. In state 2, particles undergo thermal diffusion 
along the filament. Steppers can switch from state 1 to 
state 2 with probability p\<z in a region of width e around 
the potential minimum. The time spent diffusing in state 
2 is determined by the homogeneous probability pai to 
jump back to state 1. For energy consuming particles, the 
ratio P12/P21 differs from its thermal equilibrium ratio, 
thus breaking detailed balance and leading to rectified 
motion along the filament (sec Fig.[T]) flij . 

To account efficiently for hydrodynamics that capture 
the geometrical confinement in which steppers displace, 
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FIG. 2: Average velocity of steppers as a function of particle 
filament coverage, For comparison the corresponding 

results for "sliders" (see text) are also shown: open and filled 
down-triangles. Inset: comparison between R=l (triangles) 
and R=0.8 (squares) steppers including HI. 



the fluid surrounding the filaments and moving particles 
is modeled as an ideal dissipativc fluid where point parti- 
cles closer than a prescribed cutoff distance, r c , interact 
through the Lowe- Andersen thermostat [l2|. Fluid par- 
ticles within r c from the filament surface or moving par- 
ticles exchange momentum using the same local thermo- 
stat. The model thus enforces local equilibrium and lin- 
ear and angular momentum conservation, yet it does not 
change the free-energy landscape of the ratchet model. 
The spatial extent of moving particles is taken into ac- 
count modeling them as hard, spherical particles with 
radius a. We have studied the collective motion of step- 
pers that move at a fixed distance r from the surface of 
a straight, fixed, filament of cross-section 7ttq. The fila- 
ment is aligned parallel to one of the edges of our simula- 
tion box, of length L ~ 100 1 and square section with edge 
L x = L y = 24 1 and periodic boundary conditions are 
applied (Similar results have been obtained for cylindri- 
cal confinement). For simplicity, moving particles never 
detach from the filament (i.e. they are infinitely proces- 
sive according to the classification of molecular motors). 
Starting from a particle random distribution along the 
filament, we monitor the stepper average velocity v, and 
its dispersion on time scales during which moving parti- 
cles can cross the periodic unit box (t < L/vq), being v 
the average velocity of an isolated particle. 

Before turning to hydrodynamic interactions (HI), we 
first briefly consider the role of excluded-volume inter- 
actions (EV) among steppers. The effect of EV on v 
depends on the ratio R = a/ 1. For biological molecular 
motors R is of order one, but in artificial ratchets it can 
be tuned by controlling the range of particle-particle in- 
teractions [13j. Figure [2] shows the variation of v with the 
moving particle filament coverage. The results in Fig. [2] 
were obtained when particles move along a single track 
(prescribed by a confining angular potential) but this as- 



sumption is not essential: similar results are obtained if 
particles are allowed to spread on the filament surface. 
The mean velocity is sensitive to the precise value of the 
ratio R. A monotonic increase of v with filament cov- 
erage (4>id = 2aN/L 1 with N the number of particles 
on the filament) is observed for incommensurate moving 
particles, opposed to what is observed for passive diffu- 
sion. This behavior can be understood because when the 
range of particle-particle repulsion is not commensurate 
with the ratchet period, a stepper in the diffusive state 
may be pushed toward the ratchet maximum (i.e. to the 
right into 5 in the inset of Fig. []} by a neighbor stay- 
ing on its left (in region 7). Hence, steppers speed up 
due to the decrease in the time it takes them to move 
to the next ratchet minimum. When the particle size 
is commensurate with the ratchet period, R = 1, v de- 
pends only weakly on filament coverage, except at very 
high concentrations where many-particle ratchet motion 
becomes less efficient than isolated particle motion. The 
motion of tightly elastically coupled motors also displays 
an analogous dependence on the commensurability of mo- 
tor separation with the ratchet period 141 ? ]. 

To focus on the effect of HI, we consider independent, 
processive steppers of size R = 1, where EV are unim- 
portant. Fig. [2] shows the dependence of v for particles 
walking along a filament, as a function of the fractional 
filament particle coverage, 4>\d- For (put < 0.1, the av- 
erage particle velocity does not depend on concentration 
but at higher concentrations, we observe a fairly linear 
dependence of v on 4>id- The velocity increase is signif- 
icantly larger than that observed for hydrodynamically 
coupled particles that move on a surface under the in- 
fluence of a constant external force that has been cho- 
sen such that it reproduces the average velocity of iso- 
lated particles (referred to as "sliders" in Fig. [2]) Q . The 
key difference between sliders and the ratchet motion 
characteristic of steppers is that sliders move smoothly, 
rather than in bursts. A simple mean-field argument 
can be used to estimate the particle density at which 
hydrodynamic speed-up becomes significant. The aver- 
age fraction pp± of bound particles which move under 
the influence of the ratchet force, /, (state 1) induces 
a mean drift velocity v ~ ^e^f^PPl^ j^R 2 on tnc 
diffusing steppers (state 2) over the characteristic time 
At ~ ^ /(T-<5) * n w hich bound particles displace along 
the filament. Hydrodynamic correlated motion plays a 
role when the mean displacement felt by the diffusing 
particles due to the ratcheted ones allows them to sur- 
mount the characteristic potential barrier, i.e. vAt > 6. 
Hence, hydrodynamic speed-up becomes significant for 
> 4>ie> = 28/(1 — d) 2 InL/Aa, which corresponds to 
~ 0.1 for 5 = 0.2, R = 1, L = 100 and Pl ~ 1. 

This hydrodynamic coupling, independent on the di- 
rect forces steppers exert on each other, is qualitatively 
different from the collective motion of particles interact- 
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FIG. 3: Left: Correlated (uncorrelated) motion of three ad- 
jacent particles with HI (EV). Right: probability distribution 
of the dwell time, St, between the consecutive stepping of 
neighbor steppers. 
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FIG. 4: Hi-induced clustering of steppers moving on the cylin- 
drical surface of the filament. Filled symbols show the aver- 
age velocity of steppers (triangles) and fluid (circles) before 
clustering. The corresponding open symbols show the veloc- 
ity after clustering. "Sliders" (downward triangles) clustering 
does not affect the mean velocity. 



ing through short range forces where the possibility to 
enhance collectively their velocity depends on the degree 
of attraction and its commcnsurability with the poten- 
tial ratchet 0, EH- When (f)m approaches one, EV lead 
to a substantial decrease of the mean particle velocity. 
Fig. [5] also shows that collective motion induces a net 
average flow of the fluid in which the filament and parti- 
cles are embedded. Hence, active ratchet motion favors 
also the transport of suspended, passive particles. This 
transport scenario is likely to be relevant for cargo mo- 
tion in elongated geometries found in biological systems 
such as neurons and in plant cells. If the particle size is 
not commensurate with the ratchet period, R < 1, the 
average stepper velocity is further enhanced, as shown in 
the inset of Fig. [5] 

The probability distribution function (pdf) of dwell 
times, 6t, between subsequent stepping events of nearby 
particles provides insight into the mechanism of hydro- 
dynamic velocity enhancement. As shown in fig. |3l when 
hydrodynamic coupling is taken into account a sharp 
peak for small 5t appears indicating that HI favors the 
correlated motion of subsequent particles. When a parti- 
cle steps, it generates a transient flow that pushes (pulls) 
the particles that are in front (behind) it, thus facili- 
tating their crossing of the ratchet barrier. A similar 
speed-up has been described for colloids moving in a saw- 
tooth potential, under the influence of a constant, exter- 
nal force [14|. Such behavior is not observed for EV in- 
teractions that, for commensurate motors (i? = 1) show 
a behavior very close to the, reference, behavior of free 
motors (nor HI nor EV interaction), while for smaller R, 
the short-ranged interactions between moving particles 
hinder their relative displacements, which favors shorter 
dwell times. The left panels of Fig. [3] display typical tra- 
jectories of three nearby steppers; in the presence of HI 
(top), jumps are strongly correlated, while no correlated 
jumps are observed in the absence of HI (bottom). 

Thermal ratchet motors are not necessarily restricted 
to move along a single track. For example, in biological 
systems microtubules indeed are composed by a num- 



ber of polarized tracks arranged in a cylindrical filament. 
To mimic the latter scenario we allowed particles to dis- 
place freely on the filament surface while subject to the 
same filament interaction, V(x). An even stronger in- 
crease in v is observed now as a function of the surface 
particle coverage, (j>2D = Na 2 /(2rL), as shown in Fig. 2) 
The more homogeneous filament coverage with respect 
to single track motion leads to average particle velocities 
that are one order of magnitude larger and it results in 
a stronger coupling to the surrounding fluid; the average 
flow velocity is comparable to that of the steppers. This 
implies that collective active motion can induce efficient 
passive intracellular transport. As in the single track 
case, there is a characteristic threshold coverage, <\>id-, 
above which hydrodynamic speed up of v is observed. 

The fact that the mean particle velocity depends on 
the stepper concentration leads to "bunching" . A group 
of moving particles with a higher than average concen- 
tration will move faster than the rest, thus leading to 
the build up of larger clusters. This clustering results in 
heterogeneous filament coverage. We have analyzed clus- 
ter formation of stepper aggregates on long time scales, 
t ^> L/vp. Using a robust distance criterion to identify 



clusters [17J, we find that at intermediate and high fil- 



ament coverage particles aggregate into a single cluster- 
that survives for the rest of the simulation up to twice the 
time needed to form the aggregate (see panel A of Fig. [5]). 
At long times the average particle speed and the coupling 
with the embedding fluid decrease slightly, as shown in 
Fig. @J underlining that the heterogeneous coverage of 
the filament jointly with an increased effective filament 
coverage caused by clustering affect the particle collec- 
tive motion. As the mean filament coverage increases, 
the aggregate percolates along the filament. At low fil- 
ament coverages, 4>2D < 0.15, we do not observe large 
aggregates. Rather, steppers organize into small clusters 
that form and dissolve. However, we cannot rule out that 
clusters develop at longer times. 

Clearly, the hydrodynamic speed-up of moving par- 



4 




^yvTip^ll' il 'm 







10 20 
t/T 



D =0.350 
D =0.175 
D =0.094 
D =0.044 

_J i 

30 



0.2 0.4 0.6 0.8 



10 



FIG. 5: Left: fraction of steppers belonging to the largest 
cluster, as a function of time. The time unit r = L/vo- 
Right: Average number of stepper power strokes to displace 
one ratchet period. The energy consumption is proportional 
to the number of power strokes. Up (Down) triangles stands 
for steppers interacting through EV and R=0.8 (R=l) while 
filled (open) circles stand for steppers interacting through HI 
for single track (free diffusion on the filament surface) motion. 



ticle groups increases transport efficiency. Assuming 
the transition rate between the two internal states of a 
particle is independent of filament coverage, the energy 
consumption of non-interacting steppers is constant. 
A good measure of particle displacement efficiency is 
then the average number of cycles, Nj, between the two 
internal energy states required to make a particle move 
one period. Panel B of Fig. [5] shows that Nj decreases 
sharply when the steppers speed-up due to HI. For EV 
the decrease in Nj is much more gradual and it is essen- 
tially absent if the size of the particle is commensurate 
with V(x). The sharp decrease of Nj with coverage 
illustrates that, as a result of hydrodynamic coupling, 
many particles "surf along" on the flow field generated 
by the power stroke of a single stepper. 

We have shown that dynamic interactions among 
weakly coupled steppers strongly affect their average 
speed due to the changes in the particles' dwell times 
between successive jumps when the particle motion pro- 
ceeds in short steps. Such an interaction leads to a speed- 
up of up to a factor 60 compared with the velocity of an 
isolated thermal ratchet motor. Moreover, HI induces a 
surrounding fluid flow that allows transport of suspended 
particles with a mean velocity comparable to that of par- 
ticles driving the flow. The phenomena that we describe 
should be relevant and experimentally observable, for ex- 
ample, in the collective transport of colloidal particles 
moving under the action of a ratchet potential gener- 
ated by an optical trap although to our knowledge, such 
experiments have not yet been attempted. We expect 
our results are also relevant for biological systems where 
the steppers are molecular motors pulling a cargo. In 
particular, the geometry studied is of particular interest 
for axonal transport [{§ [Io[ where the proccssive motor 



KIF1A has been proposed as a good experimental real- 
ization of the two state ratchet model. More generally, 
we expect our results to be relevant whenever the het- 
erogeneous dynamic described by the two state model is 
fulfilled. This is found , for example, of non-processive 
motors whose tails arc embedded in a membrane-covered 



cargo. In this case the fluid-like membrane [llj will sus- 
tain the hydrodynamic coupling while the proximity of 
the cargo to the filament will enhance the binding rate 
of non-processive motors. Finally we speculate that the 
strong velocity enhancement that we observe may play a 
role in the very fast cytoplasmic streaming that has been 
observed in some plant cells [|| . 
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